Lightly doped La2_ x Sr x Cu04 as a Lifshitz helimagnet 
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We study the static magnetic correlations in lightly doped I^-zSr^CuCU within the framework 
of a dipolar frustration model for a canted antiferromagnet. We show that the stability of the canted 
Neel state for x < 2% is due to the Dzyaloshinskii-Moriya and XY anisotropies. For higher doping 
the ground state is unstable towards a helicoidal magnetic phase where the transverse components of 
the staggered magnetization rotate in a plane perpendicular to the orthorhombic b-axis. Our theory 
reconciles, for the first time, the incommensurate peaks observed in elastic neutron scattering with 
Raman and magnetic susceptibility experiments in I^-iSr^CuCU. 
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PACS numbers: 74.25.Ha, 75.10.Jm, 74.72.Dn 

Introduction.— The magnetic and transport proper- 
ties of La9- a; Sr a ;Cu04 (LSCO) evolve tremendously with 
doping [y. For < x < xaf ^ 0.02, LSCO is a 
Mott-Hubbard insulator that exhibits 3D antiferromag- 
netic (AF) order with a Neel temperature Tjv — 325K at 
x = 0. For xaf < x < xsc — 0.055, inside the spin- 
glass (SG) phase, static incommensurate (IC) magnetic 
correlations are observed in neutron scattering, demon- 
strating the importance and persistence of the strong AF 
correlations well after the Neel long range order has been 
destroyed. Further doping, x > x$c, eventually drives 
the system into a superconducting (SC) state, where dy- 
namic magnetic correlations are also known to play an 
important role pj. 

Inelastic neutron scattering (INS) experiments within 
the SC phase have revealed that dynamical IC spin cor- 
relations coexist with superconductivity Q- Moreover, 
the experimentally observed four peaks associated with 
the IC magnetic order are accompanied by charge peaks, 
with periodicity twice the magnetic one, provided a low- 
temperature tetragonal phase is stabilized (as, for ex- 
ample, by Nd doping) j3j. These features have been 
consistently interpreted within a model of stripes act- 
ing as antiphase domain walls 0. More recently, static 
IC magnetic correlations have also been observed within 
the low-temperature orthorhombic (LTO) SG phase of 
LSCO 0|. In this case, only two IC peaks, rotated by 
45° with respect to the ones in the SC phase, have been 
observed. Because such incommensurability follows the 
same linear dependence upon doping as the other four 
peaks @, the above results, as well as recent Raman 
measurements Q , have been interpreted as a signature of 
diagonal stripe order. However, associated charge peaks 
have never been observed in the LTO phase of LSCO, 
what raises the question of the validity of the diagonal 
stripe picture within the insulating SG phase. 

Recently, a model based on the seminal ideas of 
Shraiman and Siggia [£| has been proposed. It ac- 
counts for the observed magnetic IC diagonal peaks with- 
out having to assume a charge order Within this 
model, the two IC magnetic peaks would originate from 



a spiral state, with the staggered magnetization rotating 
within the Cu02 layers (ab plane), in such a way that 
the measured magnetic incommensurability would corre- 
spond to the inverse spiral pitch. More recently, Sushkov 
and Kotov have considered a similar spiral state in the 
t— t' — t" — J model ^J]. It was shown that while in the in- 
sulating region, x < xsc, the diagonal (1, 1) and (1, —1) 
spirals have lower energy than the horizontal (1,0) and 
vertical (0, 1) ones due to the Coulomb trapping of the 
doped hole near the Sr ion, in the metallic state, x > x$c, 
the Fermi motion energy favors the (1, 0) and (0, 1) spiral 
states, leading to a jump of 45° in the direction of the 
spiral pitch across the metal-insulator transition. 

Despite being able to qualitatively explain neutron 
scattering data in the SG phase, the above spiral pic- 
tures have two major problems: i) the collinear Neel state 
is unstable towards the spiral state already at infinites- 
imal doping, x ^ 0; ii) the spiraling of the staggered 
magnetization in the Cu02 layers is not consistent with 
magnetic susceptibility experiments in LSCO, which indi- 
cate that the Cu ++ spins remain confined to the be plane 
throughout the SG phase (Xa remains featureless) flll |. 
In this Letter we demonstrate that only by considering 
properly both the Dzyaloshinskii-Moriya (DM) and XY 
anisotropies the above two problems can be solved, and 
we provide an elegant description of the unified physics of 
IC correlations and anisotropic magnetic response in the 
SG phase of LSCO. We show that the anisotropies give 
robustness to the canted Neel state for x < xaf and that 
the same anisotropies lead to the formation of a Lifshitz 
helimagnet above xaf, with the small transverse compo- 
nents of the staggered magnetization rotating in the ac 
plane. Such helix rotation of the staggered magnetization 
is accompanied by a small precession of the local magne- 
tization (weak-ferromagnetic (WF) moment) around the 
out-of-plane c-axis. The space integrated sublattice mag- 
netization and WF moment, on the other hand, are ori- 
ented along the b- and c-axis, respectively, in such a way 
that, from the point of view of the magnetic susceptibil- 
ity, the total spin is confined to the be plane, in agreement 
with experiments (lll |. 
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The model— The presence of static magnetic correla- 
tions within the SG phase of LSCO allows us to use a 
Hamiltonian description of the long-wavelength fluctua- 
tions of the staggered order parameter, n (see 0,E3] for 
a detailed derivation of the model). For x = we can 
write (we use a soft version of the n 2 = 1 constraint) 

(i) 

where n is the local staggered magnetization, t 1 = p s /T 
is the renormalized classical spin stiffness, with p s — JS 2 , 
J is the AF superexchange, a is the lattice constant, T is 
the temperature, and S = 1/2 is the spin. The effect of 
both the DM and XY anisotropies is to generate gaped 
transverse excitations. We define m a , (a = a, b, c), as the 
bare masses of the magnetic excitations along the three 
directions of the LTO phase. The bare masses m a and 
m c are related to the DM and XY anisotropy parame- 
ters, respectively, whereas the mass of the longitudinal 
mode, to;, (the heaviest of all three), is proportional to 
the coupling constant uq. Furthermore, the uniform part 
of the Cu ++ spins, L, is related to the staggered order 
parameter, n, through (L) = (l/2)((n) x D + ), and gives 
rise to a WF moment perpendicular to the CuC>2 layer 
[HQ. Here D+ = (D+,0,0) is a thermodynamic DM 
vector oriented along the a-axis, related to the tilting 
angle of the oxygen octahedra D + ~ 6 and we use 
units such that J = a = I. 

The effect of a small number of holes (x ^ 0) on an 
antiferromagnet can be characterized, in a coarse-grained 
description, in terms of a dipolar frustration of the AF 
background 0,0. We introduce a dipolar field, P M , that 
couples to the AF magnetization current as 
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2A / d 2 xP M • (n x <9 M n) 



(2) 



where A = X/T is the dipolar coupling constant, A ~ 1 
[l^j . and the field is a vector in both lattice and spin 
spaces. The above mathematical form belongs to a class 
known as Lifshitz invariants, which play an important 
role in stabilizing long-period spatially modulated struc- 
tures with fixed sense of rotation of the vectors n in sys- 
tems like Ba 2 CuGe 2 7 and K 2 V 3 8 0. 

In the case of a nonuniform configuration for the dipo- 
lar field, P M , we can write 

HD = hj ^{(vp^+a^op;) 2 }, (3) 

where K = kT with k denoting the dipole stiffness, and 
fi afJ , is the bare mass of the dipole field, P", which, in 
microscopic terms, is related to the energy cost of pop- 
ulating the pth valley of the vacancy Fermi surface with 
the spin polarization a 0, Q] . 

The appearance of elastic IC peaks in neutron scatter- 
ing can occur, in principle, whenever (P M ) ^ 0. We thus 



calculate the effect of the magnons on the self-energy of 
the dipolar fields, arising from the dipole-magnon inter- 
action in Eq. Q. Straightforward calculations yield 
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The Fermi wavevector kp = ^/nx provides the cutoff in 
the momentum integrals, because the coarse-grained de- 
scription employed here is valid on length scales much 
larger than the distance between the holes, I ^> kp - 
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0, Q|. Inspection of Eq. <@J reveals that, even 
when the bare dipole masses are isotropic in the spin 
space, = Va, the renormalized ones may be 

anisotropic due to the dipolar interaction, Eq. |0), which 
involves different magnon modes. Moreover, an instabil- 
ity to a phase with (P") ^ may indeed occur. In par- 
ticular, when m a = m c — 0, such instability appears al- 
ready at infinitesimal doping, because the momentum in- 
tegral in Eq. <@J diverges logarithmically. Finally, due to 
the hierarchy between the masses of the magnon modes, 
m a < m c <C nib (actually m& — > oo in the nonlinear 
n 2 = 1 description), we find that M 2 ^^ will become nega- 
tive first for the spin component of the dipolar field along 
the orthorhombic 6-axis. 

Canted Neel state.— Let us consider first the low dop- 
ing regime where (P™) = 0. In this case, the magnetism 
is commensurate and the fluctuating dipoles give rise to 
corrections to the magnon self-energy in one-loop ap- 
proximation. The renormalized propagator for the n a 
field reads G a| a(q) = tS a f3(q 2 + 77i 2 /2 — E a (q)) -1 , where 

E a (q) = K(txy 

d 2 k {k + q)n(k + q)n 



(2tt) 2 (fc 2 + m 2 /2)[(k-q) 2 + M y 

The above self-energy is an analytical function of the 
momentum and thus we write E a (q) = E a (0) + 
(l/2)q tx q l/ d 2 ^ a (q)/dq fl dq^\q =0 . The first term in the 
self-energy, E Q (0), gives rise to corrections to the magnon 
mass 
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whereas the second term leads to a reduction of the spin 
stiffness, 



d 2 k {k 2 ) 2 +^ v 



(2tt) 2 (k 2 + m 2 /2)(k 2 +p 2 )3 



(6) 



It is clear now from Eqs. 10} and © that it is exactly 
the DM and XY anisotropies that give robustness to the 
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canted Neel state at low doping, because they lead to 
nonvanishing magnon gaps, m a ^ and m c ^ 0. Most 
remarkably, we find that the reduction of the DM gap, 
M a , is consistent with recent Raman scattering experi- 
ments by Gozar et al. ^f|, which show that the DM gap 
decreases with doping and vanishes at x = xaf — 2%, 
when the SG phase sets in. Experimentally, at x = 1% 
the decrease of the DM gap is about 26% [lfl |. Such 
huge renormalization, however, cannot be explained sim- 
ply from the decrease of the tilting angle of the oxygen 
octahedra in the LTO phase of LSCO, which is known 
to happen with Sr doping. Our model yields the de- 
crease of the DM gap with doping due to the coupling of 
the background magnetization current to the dipoles, Eq. 
<EJl- Using Eq. (0 with the bare values of parameters, 
k ~ 1.1, A ~L, t ~ 1/S 2 , m a ~ 2.5 • 1(T 2 , m c ~ 5 • 1(T 2 , 
and ji ~ 1 |l4l |. we determine the doping dependence of 
the DM gap, at T = 10 K (see Fig. EJ. The reduction 
of the DM gap at x = 1% is about 27%, in remark- 
able agreement with the experiments. Moreover, from 
the vanishing of the DM gap we obtained the critical 
concentration x^p — 0.0198 ~ 2%, as experimentally 
observed signaling the instability of the Neel phase. 




Sr concentration x (%) 

FIG. 1: Doping dependence of the DM gap as given by Eq. 
10 using the bare parameters of the model. Experimental 
data from Ref. taken at T = 10 K. 

Our results are also consistent with recent magnetic 
susceptibility measurements Within our model, the 
only effect of doping is the reduction of the magnon 
masses and of the spin stiffness (see Eqs. (0 and JSJ )• 
Therefore, we expect the qualitative features of the sus- 
ceptibility in the AF region of LSCO to remain the same 
as in the undoped compound The peaks in both 

Xb and Xc susceptibilities will be shifted towards lower 
temperatures, because both the reduction of the magnon 
gaps and of the stiffness with doping lead to a decrease 
of the Neel temperature, Tn- This is indeed observed ex- 
perimentally [llj. Furthermore, we find that the unusual 
T = hierarchy of the susceptibilities, Xa < Xc < Xb, 



is preserved |l3| (once the anisotropic van Vleck contri- 
bution to Xc is subtracted with Xa ~ &o/t, Xb = 
{l/t){D\/M*{x)), and Xc « Xa + (lAX^M^)), 
where cr measures the Cu ++ effective moment 




FIG. 2: (Color online) Helix rotation of the staggered magne- 
tization for x > xaf- Notice that the longitudinal component 
is the largest, and the small transverse ones describe an ellipse 
in the ac plane. 

Helicoidal state.— Let us consider now the case of 
higher doping where an instability towards a phase with 
(P") 7^ may occur. The easy axis for the spin part of 
the dipolar field is the orthorhombic 6-axis, which is a 
consequence of the hierarchy m a < m c <C rrib. As the 
lattice part is concerned, our model is isotropic, and, 
as such, cannot determine the lattice direction of the 
dipolar field. Motivated by neutron scattering exper- 
iments, we assume that the dipolar field acquires the 
expectation value along the 6-axis in the lattice space. 
In order to obtain the corresponding nonuniform config- 
uration for the staggered order parameter we consider 
(P M ) = (0, Po, 0)6^, and solve the equations of motion 
for the n fields. The resulting local staggered magneti- 
zation reads 

n(x) = {a a cos (Q • x), a b , a c sin (Q • x)), (7) 
with a c = -cr a (Q 2 + m 2 )/P (3 and 

2Q 2 = P 2 - m 2 - ml + J (Pq - m 2 - m 2 ) 2 - 4m 2 m 2 , 

(8) 

for Pq > m 2 + to 2 , where Pq = — p s XPo- The equa- 
tion of motion for the longitudinal component yields 
n 2 = 1 — n 2 — n 2 , and that determines Ob- The static con- 
figuration for the staggered order parameter is nonuni- 
form and its transverse components, n a and n c , rotate in 
a plane perpendicular to the b direction. Therefore the 
staggered order parameter becomes helicoidal, as shown 
in Fig- EI Such a magnetic texture arises as a result of the 
competition between the anisotropies, m a and m C) which 
favor the staggered magnetization along the 6-axis, and 
(P^*), which supports a rotation of the n field in the ac 
plane. Since (Pj?) ~ x is small, the resulting spin config- 
uration has small components in the ac plane, as shown 
in Fig. 
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The local magnetization will precess around the c-axis, 
see Fig. El because (L(x)) = (1/2) ((n(x)) x D + ), but 
the total space integrated magnetization, which is mea- 
sured by magnetic susceptibility measurements, is still 
along the c-axis 



M 



V 



d 2 x(L(x)> = (0,0, 



er b L> 4 



with a 2 = 1 — <J 2 /2 — <j 2 /2. Since the magnetic suscep- 
tibility is obtained from the total magnetization, M, we 
conclude that the nonuniform configuration in Eq. J7J 
is consistent with the experiments of Lavrov et al. 
which indicate that the Cu ++ spins remain confined to 
the orthorhombic be plane throughout the SG phase, thus 
yielding a featureless \a EH- 

The new helicoidal magnetic structure gives rise to IC 
peaks in neutron scattering at wave vector Q = (0, Q, 0), 
corresponding to the inverse helix pitch. Assuming that 
in the SG phase Eq. © holds with renormalized magnon 



we obtain Q 2 
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gaps, M a and M c 
because M a (x > xaf) = 0. This expression allows us 
to understand both the linear dependence of the incom- 
mensurability at higher doping x, when M c is negligible, 
as well as a possible deviation from linearity in a region 
of sizable M c (x — > xaf)- In fact, neutron scattering 
experiments do exhibit such deviation from linearity at 
x = 0.024 @. 

We can now make an important new theoretical predic- 
tion about the behavior of the IC peaks in the presence of 
an applied magnetic field B _L to the CuC-2 planes. Since 
in this case the XY gap softens [l^ j. the incommensu- 
rability is expected to reapproach the linear behavior, 
Q ~ x, and is thus larger than for the zero field case. In 
the stripe picture, on the other hand, the peak positions 
remain unchanged ^jj and thus a perpendicular mag- 
netic field can be used in order to decide between these 
two scenarios. 

Finally, we provide an independent estimate for the 
critical doping concentration x^p at which such heli- 
coidal instability takes place, from the behavior of the 
dipolar mass Mbb at T = 0. Assuming that Mbb van- 
ishes at x^ F we find (recall that kp — kp(x)) 



2 

^bb 




Mc 



(9) 



In order to estimate x^p we use the bare value for the 
mass of the XY mode at the critical point, m c ~ 5 • 10~ 2 , 
the spin-wave velocity c = 2SV2, and for lightly doped 
LSCO (a; ~ 0.01), we have ck F /M c ~ 5, with kp = 
y/irx. Thus, in Eq. Q we can consider M c jckp as a 
small parameter to obtain x^p 



0.0206 ~ 2.1%, since it ~ 1.1, A ~ 1, p s ~ S 2 , and 
Mbb ~ 1 01- consistent with the vanishing of the DM 
gap, and in agreement with experiments. 

Conclusions.— We propose a description of lightly 
doped LSCO in terms of a dipolar frustration model of a 
canted antiferromagnet. For x < xaf the dipole-magnon 
interaction leads to a reduction of the DM and XY gaps, 
as well as of the spin stiffness, and the qualitative fea- 
tures of the magnetic susceptibility remain the same as 
in the undoped compound |l3|, in agreement with exper- 
iments [Ti!|. The robustness of the AF state stems from 
the DM and XY anisotropies, and for xaf < x < xsc 
the ground state is unstable towards a helicoidal phase 
where the local WF moment precesses around the c-axis 
(see Fig. |2j. Such helicoidal magnetic structure gives 
rise to two IC peaks along the b direction, as observed 
in neutron scattering while yielding a featureless \ a 
fllj | . The incommensurability is expected to scale linearly 
with doping, but to deviate from linearity as x — > x\ F . 
We also predict that such linearity is recovered once a 
perpendicular magnetic field is applied. Our theory pro- 
vides, for the first time, a consistent description which 
reconciles neutron scattering, Raman, and magnetic sus- 
ceptibility measurements in lightly doped LSCO, and in- 
serts the latter into a wider class of Lifshitz helimagnets 
such as Ba2CuGe207 and K2V3O8 [lfij - 
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